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A T-FUNCTION SOLUTION 
TO THE SIXTH PAINLEVE TRANSCENDENT 

YU. V. BREZHNEV 

Abstract. We represent and analyze the general solution of the sixth Painleve tran- 
scendent in the Picard-Hitchin-Okamoto class in the Painleve form as the logarithmic 
derivative of the ratio of certain T-functions. These functions are expressible explicitly 
in terms of the elliptic Legendre integrals and Jacobi ^-functions, for which we write the 
general differentiation rules. We also establish a relation between the V& equation and 
the uniformization of algebraic curves and present examples. 
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1. Introduction 
The first appearance of the sixth Painleve transcendent 

n : y xx = ^[- + — 7 + )y x - - + + )y x 

2 \y y — l y — xj \x x — l y — x/ 

, y(y - - x) ( x x-i ( i^x(x-i) 



( a x , x ~ 1 Ik l \ 

( v a - /3 ^ +7 (^T)2-l S -2) 



x 2 {x — l) 2 \ y 2 (y — l) 2 ' V'ijj — x)* 

dates back to the 1905 Fuchs' paper [9]. Contemporary applications of the equation 
are very well known [3, 25, 13] and the most nontrivial ones are in cosmology. Fuchs had 
already used the elliptic integral in his paper, and Painleve a year later wrote this equation 
in the remarkable form 

7T 2 d 2 Z 

-T^2= ap'(z\T) + (3p'{z - l|r) + jp'(z - r\r) + 5p'{z - 1 - r|r), (2) 

again using elliptic functions and performing the transcendental variable change (y, x) i— > 
(z,r): 

<%tY V 3 + 3^(r) vr 2 ^|(r) W 
(these formulae are equivalent to the original Painleve ones; see relations (4) below). Here 
and hereafter, p(z\r) := p(z\l, r) and other Weierstrass' functions (a,C,P,p') [lj 8, 26] 
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are constructed with the half-periods (u,u)') = (1, r), and the Jacobi theta constants f?(r) 
depend on the modulus r belonging to the upper half-plane H + , i.e. S(r) > (see the 
appendix for definitions). 

This Painleve result, known as the p-form of the "Pe-equation was rediscovered and 
developed in the 1990s with other approaches, and the deep relation between (1) and 
second-order linear differential equations of the Fuchsian class is well known. The tran- 
sition from these equations to Painleve representation (2) was recently well described 
informally [2]. In the same place additional references can also be found. 

The interest in general Painleve equation (1) increased dramatically after Tod demon- 
strated in 1994 [ ], that imposing the SU(2)-invariance condition on the metric (the 
famous cosmological Bianchi-IX model) together with assuming the metric conformal in- 
variance are consistent with the Einstein- Weyl self-duality equations R a b = ^Ag^ and 
Wj bcd = 0. These equations are partial differential equations but they can be reduced 
to an ordinary differential equation with respect to the conformal factor, which can be 
expressed in terms of a solution of the V% equation for particular values of the parameters 
(q, /3, 7, 5) (sect. 2.1). Equation (1) with arbitrary parameters was later obtained as a self- 
similar reduction of the equations for the Ernst potential in general relativity [23]. This 
is not the only example of the appearance of equation (1) in applications: it suffices to 
recall the Yang-Mills equations (Mason- Woodhouse (1996)), two-dimensional topological 
field theories (Dubrovin), etc. In the general setting equation (1) was also derived in [17] 
as a self-similar reduction. See also references in these works. 

1.1. The Painleve substitution (3). Equation (1) and substitution (3) reflect the prop- 
erties of not only the equation itself but also its solutions. If an equation is of the 2nd 
order, linear in y xx , and rational in y x and if its solutions have only fixed branching points 
(the Painleve property), then the number of these points, as is known, is reducible to 
three. We can always choose them to be located at the points Xj = {0,1, oo}, and the 
equation then either takes Fuchs-Painleve form (1) or represents some of its limit cases 
[20]. In turn, the plane of the variable x can be conformally mapped one-to-one to the 
fundamental quadrangle of the modular group T(2) [8] in the plane of the new variable 
r using the modular function x = A/ 2 (r), i.e. the first formula in (3). Function k l2 {r) 
behaves exponentially with respect to the local parameter r at the preimages of the points 
Xj. Because of this the branching y = y(x) of an arbitrary power-law or logarithmic char- 
acter in the vicinity of Xj transforms into a locally single- valued dependence on r [4, §7]. 
This partially explains 1 the origin of substitution (3). Fuchs obtained this substitution [ 
using the Legendre equation and Painleve used the modular function [20]. 

Painleve himself wrote the equation not in form (2) but in terms of the modular function x = with 
the three singular points Xj = {0, l,oo}, appearing when inverting the elliptic integral 



This is a partial explanation because a general rigorous statement on the branching of solutions of the 
Va equation has apparently not yet been introduced in the literature. For example, it is unknown whether 
it is (at most) a local branching of the type of the rational function x a \n n x. See the important paper [12], 
which was devoted to precisely this question and where Fuchsian equation (5) — an equivalent to equation 
(2) — was the main object of study. 
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* = l[ V /( tv V » = '<^ (4) 

w i Joo\/y{y -l){y-x) ^ 

In his original notation, the equation Va in the p-form is as follows [20, p. 1117]: 

where the misprint 77 h> z should be corrected. Equation (1), has the form of an 'inhomogeneous' Legendre 
differential equation [9] 



d 2 u 2t - 1 du u _ \J A(A — 1)(A — t) 

dt? + t(t -i)dt + 4t(t - 1) ~~ 2t 2 (t-l) 2 



i f - 1 t(t- 

"Joa &o To H - "'l 7^ ™t 



A 2 (A — l) 2 '(A-t) 2 



(5) 



where it = f . dx and A = X(t) is a solution of the equation Vr. 

1.2. The motivation for the work and the paper content. Other singular points are 
movable but are poles, and the global behavior of solutions is therefore totally governed 
by the logarithmic derivatives of 'entire' functions [20], [7, pp. 77-180]. Painleve stressed 
that construction of such transcendental functions completes the procedure of integrating 
an equation (the concept of the 'integration parfaite' in the Painleve terminology). The 
pole character is known for all Painleve equations [11, 7] and all their solutions except 
those for the equation V\ have the structure [20, p. 123] 

y ^ L n _2 (6) 

y dx 7i w 

(in the generic case) with the functions T 12 (x) having no movable singularities. See the 
works [ , p. 371] and [7, p. 165] for the complete list of such formulae. 

If we disregard the known automorphisms of solutions of the equation in the space of 
the parameters (a, /?, 7, 5) (see, e.g., the Gromak paper in [7], [15], [11, §42], [14] and the 
references therein), only two cases are known where the general integral for equation (1) 
can be written. These are the solutions of Picard [22, pp. 298-300] and Hitchin (sect. 2.1), 
but Painleve form (6) is known for neither of these solutions. 

On the other hand, representing solutions in terms of the corresponding analogues of 
T -functions occurs in almost all investigations of nonlinear equations somehow related 
to the integrability property: in the Hirota method, in the theory of isomonodromic 
deformations [15], in soliton theory, and in its O-functional generalizations. In this respect, 
it is interesting to show that the Painleve equations also admit a nontrivial situation in 
which the functions T x 2 { x ) can be written explicitly when we have the general integral of 
motion. Here, we fill this gap (sects. 2.2, 4) and also describe the distribution of solution 
singularities (sect. 3). Technically, we need special differential properties of Jacobi's theta- 
functions and their theta-constants. These properties are also new, and we present them in 
the appendix. As follows from the solution (this was properly noted in [15]), the character 
of the Painleve T -functions differs drastically from that of the O-functional solutions of 
the solitonic equations. Nevertheless, a variant of the straight-line section of a manifold is 
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present (sect. 4). In section 5 we demonstrate the connection of the equation V% with the 
uniformization of nontrivial algebraic curves. To what was said in sect 1.1, we add that 
Painleve himself explicitly wrote [7, p. 81] about representing solutions in terms of single- 
valued functions in the context of the Painleve property. Picard mentioned the same thing 
many times in [21, pp. 93, 188], [22]. The relation to uniformization mentioned below and 
the presentation of numerous examples (sect. 5.2) is therefore not accidental. 

2. The Painleve form 

2.1. The Hitchin solution. The Hitchin solution corresponds to the parameter values 

a = (3 = 1 = 5 = ^, (7) 

which were found in [13] under the condition that the Einstein equation metric is conformal 
with SU(2)-invariant anti-self-dual Weyl tensor [25]. Equation (2) in this case can be easily 
reduced to the equivalent form written in theta functions 

where r] is the Dedekind function 

oo 

rj( T ) =e^ T JJ(l-e 2 ^ ir ). 

fc=i 

The general solution to equation (8) in parametric form (3) is [13] 

/ , n , . ,x 1 P'(At + B\t) , , 

p(Ar) = p(Ar + B\t) + - ^ - |1 _ ^L> _ ^ , (9) 

where A, B are the integration constants. This solution in a different implicit form was 
mentioned by Okamoto [19, p. 366]. 

Solution (9) is not only the most nontrivial of all those currently known but is also rather 
instructive for P-equations in general because it provides the general integral. It has been 
considered in detail in the literature and is certainly more than 'just a solution'. The Picard 
solution is degenerate in this respect because it is a perfect square (see sect. 4.1). This 
explains why we choose the Hitchin solution for further analysis. Other particular cases are 
known that contain variations of the Painleve c^Ln-form but do not contain T -functions. 
Allowing more freedom in choosing the parameters (a, f3, 7, 5), these solutions nevertheless 
contain only one integration constant and are related to the linear hypergeometric 2F1- 
equation [ , p. 733], [ I, §44], [19, p. 374 §5.4], [ , p. 145]. 

2.2. The Painleve form. One of the theta-function versions of solution (9) was proposed 
in [13] and another form in [15]. But the complexity of these solution forms was mentioned 
in [13, p. 75] itself, in [15, p. 901], and also in work [3], consisting entirely of calculations. 
In those papers, solutions were also presented in the parametric forms and contained the 
set of functions 8, 9', 8", 6"', 1?, . However one can obtain the solution in a compact 
form if we transform Weierstrass functions in (9) into the theta functions (cf. [13, p. 74]). 
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Proposition 1. The general solution to equation (1) with parameters (7) is 



9f X&x + TTiAe- 1 



(10) 



71 

w/iere 9 = 0(\At + |.B|t) and $ = i?(r). 

The possibility of representing solution (10) in explicit form (6) is related to nontrivial 
properties of the theta functions, for instance, to the identity that follows from Lemma 3 
in the appendix 

d x r >/ I X /NT /IN 1 P'( Z I T ) / \ 

2i d^^^ " ^ (T) > = ^ + 2 C(^|r) - 4r) + ^ ' 

Comparing this property with (9), we obtain the total logarithmic derivative, 

vr d^ C^T + ijlr)-^^)-^) 
Pv z r ) = ~ -j-Ln -— . 11) 

2l (XT ^ ( r ) 

Replacing (-A,-B) i— > 2(^4,-8) and transforming the right-hand side of (11) into ^-functions 
(lemma 1 sect. 7) we can write the solution in the form 

2i 1 d 9[(Ar + £|t) + 2k\A6i{At + £|t) ^(r) 
y ~7r^( r ) dr Ln 4*( T )6iAT + B\t) ' ^"^(r)' 1 j 

To rewrite the solution in the initial variables (x,y), we must use the formulae inverting 
substitution (3) and its differential. They can be written in terms of the complete elliptic 
integrals K, K' [1, 27]. It is easy to see that this transition together with the differentiation 
is determined by the formulae 

j- = nix(x - l)0&r)-f , *!(r) = - VT^K'(Vx), (13) 

UT dX IT 



and inversion of the first formula in (3) is given by the expression (series) 

6 * = ^{-4^}= 



2 



+ 8 ^ +84 —— + 992 —— +12514 —— + 



1 



16 / V 16 / V 16 / V 16 7 V 16 

Simplifying the obtained expressions and replacing \A \— > A, we obtain the sought solution 
form. 

Theorem. The general integral of equation (1) with parameters (7) has the form 

v = 2,(1 - x) * Ln m±m ± ^* ±m . (15) 

Here and occasionally hereafter, we use the shorthand notation K := K{yfx) and K' := K'(y/x) . 
The quantities K, K' can be written in terms of the hypergeometric functions [1, 8] 

K(y/x) = \-iF x (i i; l|x) , K\Jx) = \ ■ 2 F X (i, \; l|l 

but we can use these series simultaneously only in the intersection of the convergence 
domains. In applied problems, we therefore avoid using numerous transformations for 
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analytic continuations of the 2-^i-series by regarding the functions K, K' as the complete 
elliptic integrals 

1 1 

J A 2 )(l-xA 2 ) 7 y^(A 2 - x)(l - A 2 ) 

(which are everywhere finite functions) or as the Legendre functions P_i, Q_i[27]. 

Verifying this by substituting (15) directly in (1) is a rather nontrivial exercise in which 
we must use Lemma 1 and its simple corollaries. We have to use also the known dif- 
ferentiation rules for complete elliptic integrals [ ]. Analogously setting E := E(y/x) and 
E' := E'(y/x), we can write these rules as 



dx x(l — x) x dx x(x — 1) x — j 

d E K d , E' K' 

2—E = , 2— E' 



(16) 



dx xx dx x — 1 x 

and they imply the Legendre relation 



7T. „„, „,„ „„, 7T 

rj = tt] — 



' - -i & EK' + E'K - KK' 



2 2 
Speaking more rigorously, in (15), we must cancel the common 'non-entire' factor 
exp (j7rir) in the series for 6\ and 6[. But this factor and also the nonmeromorphic 
singularity of the solution determined by the factor K'(y/x) are 'fixed'. Relations (15) and 
(16) demonstrate that the solution contains an additive term, which has a fixed critical 
singularity of the logarithmic type. 

Corollary 1. The function y(x) and the solution z{t) of equation (8) have the forms 

y = ^ + 2x(l-x) ^Ln| 6 ±{A§ + B\%)+ 2tta} , (17) 
where the elliptic integral p' 1 here and hereafter is defined as 

u 

-if 1 x r dx 



To follow the meromorphic x-component (the movable poles) more closely, we substitute 
expansions of type (14) in formulae (15), (17) and in the theta-function series in the 
appendix. Weierstrass derived analogous series in relation to the solution of the modular 
inversion problem for the function k 2 (r) [26, pp. 53-4, 56, 58], [27, p. 367]. Explicit 
expansions in the vicinities of all other points can be written using the formulae in the 
appendix. Formulae (12), (15), (17) themselves provide an analytic answer for the power- 
law expansions in [13] on pages 89-92 and 108-9. 
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3. The pole distribution 
One of the two series of poles x 



of solution (15) admits an explicit parameterization, 
- ■ ./■„,„ - -.!( - : ). ra.meZ. (19) 



#4 

4 



We must supplement it with the natural restriction on (n,m) of the form Ss^ ™~^ j > 
0, which imposes no restrictions on A and B other than these numbers cannot be real 
simultaneously. Although this pole distribution is described by the simple formula (19), 
its actual form is rather involved (see Fig. 1). 





Figure 1. The poles x mn of form (19) for solutions (15), (22) under A — 
125.45 - 103.29L B = 36.710 - 69.980i and (n, m) = -30 . . . 70. 

The general and characteristic behavior of the 'orthogonal' directrices of the parame- 
terizations n = const and m = const, i.e. the pattern in Fig. 1 as a whole is shown in 
Fig. 2. The pattern in this figure can be interpreted as an analogue, although distant, of 
the two straight-line orbits of poles of any double periodic function. The pole distributions 
are deformed when the initial data (A, B) are changed, which generates various patterns, 
but we observe that all the poles concentrate near the fixed singularities Xj = {0, l,oo}, 
as should be the case. We thus obtain the character of this deformation. If we take one 
pole p — x mn under fixed (?n, ?i), then we find that it moves without changes along the 
straight-line characteristic (n + A)t — (m — B) = in the space of the initial data (A, B). 

The dynamics p = pit) for any pole of series (19) is given by the same function p(t) = -*f(i) 

but with times t k that differ by a fractional linear transformation for each pole p k = p{t k ). 
The function pit) and the structure of the fundamental domain of its automorphism group 
are well studied. Real or purely imaginary poles can also be easily described. 
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Figure 2. Two pairs of directrices n = (12, 13) and m = (129, 130) of pole 
distributions (19) for A = -195.45 + 103.29i, B = -6.710 + 79.98i. All the 
seemingly 'casual' points are also 'good' poles belonging to these directries. 



We further note that the Okamoto transformation [11, §47], [15] 

y(y - i)(y- x) 



PlCARDy !->■ HlTCHINy : %j = y + 



HlTCHIN^ I — y PlCARDy : V = V 



x(x- l)y x -y 2 + y ' 

< iv ^ (20) 
y(y - i)(y - x) 

x(x- l)y x + \y 2 - xy + \x ' 

which is known to transform the Picard and Hitchin solutions one into the other, leaves 
poles (19) invariant but creates or, under the inverse transformation, annihilates the second 
series of poles determined by roots of the equation 

((At + B\t) = A((t\t) + BC(1|t). 

Solving this equation is equivalent to finding the A-points of the transcendental function 



1 


9[ 


( A K(VS) 


- B 




2tt 


h 


( A K(^c) 


- B 





f( X ;A,B) = — ) yv ' F^f, (21) 



i.e. to solving the equation f(x;A,B) = A. The pole series determined by relation (21) 
has a much more involved description and evolution than for pole series (19). 



4. The T-functions 



4.1. The Picard solution. Pole series (19) also exhausts all the poles of solutions to 
equation (1) in the Picard case a = (3 = r y = 5 = 0. It then follows from (2) that 
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z = =>■ z = ar + b. Writing relations (3) in terms of the theta function, we obtain 
another form of substitution (3) and the equivalent form of the Picard solution in [22]: 

v __m e im ^ „ __ r jM^ki B ^Mi (22) 



A comment is in order. Strictly speaking, the original Picard solution u pic [, 2, p. 299] and the Fuchs- 
Painleve representation (1) are not identical but are related by an inverse transformation. Picard con- 
structed an example of the equation using the Legendre form Ax — \J '(1 — x 2 )(l — k 2 x' 2 ), while Painleve 
used Riemann form (4). The transition between solutions of these equations is therefore given by the 
formula 

-P 1 c(^)= 2 ^ = v / ^- ( 23 ) 
In other words, we have the property that the case of parameters a = /3 — "/ = 8 — is a case (not unique) 
of Painleve equation (1) such that the square root of any of its solution satisfies, in turn, an equation that 
also has the Painleve property, namely, the Picard equation [22, p. 299] 



d u / du\ 2 u(2xu — 1 — ail du 



/du\ 2 u{2xu 
\dx) (1 - u 2 



+ 



dx 2 \dx/ (1 — u 2 )(l — xu 2 ) dx 



1 

+ 



U(l - u 2 ) 



x(l - x)(l 



= 0, (24) 



x)(l — XU 2 ) X 

where we must amend the last term with the omitted factor |. This equation appears at the end of the 
almost two-hundred-page-long Picard treatise [22, p. 298] as equation differentielle curieuse satisfied by the 
function u — sn(aui + bcu'), regarded as a function of the elliptic modulus k 2 (= x). Therefore, the ansatz 
for studying the branching of solutions to (1) assumed in [12] 

1 + x 

y(x) = p^uj^x) + v 2 u 2 (x) + v(x; v 1 ,u 2 );u} 1 (x),u} 2 (x)) H — , 

contains a somewhat 'redundant' perfect square and can be conveniently replaced with a structure of 
Hitchin's type (15), (17), which simplifies the analysis. 

Remark 1. In addition to the Picard solution, all its transformations obtained by actions 
of discrete symmetries of equation (1) [] L, §42], for instance, by the variable change 

~_ _ ^l(r) 91{At + B\t) 

V ^ y ~Y^l~^{j)el(AT + B\T)' 

are perfect squares. Adding the modular substitutions from the group T(l) and the 
anharmonic group 53, transforming the variable x [I , §23], we obtain the complete set of 
transformations in [11, pp. 214-5] and changes of the parameters (a, /3, 7, 5). Together with 
the Okamoto transformations, all the transformations obtained above and below expressed 
in terms of the theta functions and 8, T -functions reduce to permutations of indices of the 
functions 9, 9' in formulae (10), (15), (17), and (18). By the Picard-Hitchin solution, 
we therefore understand the whole class of these solutions, without especially mentioning 
this. The Picard-Hitchin solutions are functions of \fx. On the other hand, the quantity 
y/x = s itself is related to the second-order linear differential Fuchs equation via the 
integrals K,K'{^fx) and substitution (3). This equation has four singularities located 
at the points s = {0,1,-1, 00} because the integrals K and K'(x) satisfy the Legendre 
equation with the three singularities a; - = {0, l,oo}. We thus obtain the particular case 
of the Heun equation 

1 (* 2 + l) 2 y . ^,£.0) 

2s 2 (s 2 -l) 2 Yi K' 



Yss = -~^^> Y => ^ = ^=- (25) 
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The character of quadrature (15) is therefore described by the following statement. 

Corollary 2. Painleve equation (1) in the Picard and Hitchin cases is integrable in 
(9' , 6) -function quadratures over the differential field determined by the Heun functions of 
form (25) or, equivalently, by the hypergeometric Legendre functions. 

4.2. One-parameter solutions. We note that neither formula (9) nor the first formula 
in (22) imply in any obvious way that the second-order poles from the Picard series become 
the first-order poles for (17). Analogously comparing formulae (3), (9), and (22), we find 
that the first formula in (22) does not provide the solution behavior in the limit cases for 
the parameters (A,B), while this behavior follows automatically from formula (17). For 
instance, we set At + B — > with B/A = a in (17) and (18). Taking (14) and (16) into 
account, we then immediately obtain the a-parameter family of solutions 

aE'-E + K J-k\ d T , 2 \ 

V = aK' + K ^ Z = V {(T + ah (T)}r J 

which can be usefully compared with the corresponding Hitchin formula in [13, p. 78]. 

Another case looks more complicated, and Hitchin considered it only in the expansion 
form because it results in a singular conformal structure. This case is more conveniently 
analyzed in form (11), where it corresponds to the limits At + B — > {l,r, r + l}. Setting 
At + B = et + (1 + ea) and letting e tend to zero, we obtain the solution for which we 
preserve the T -form (6): 

, . , tt d T 2(r + a)i9 2 + #2 

p (z r = — — Ln s-r 

5 v 1 ; 2 idT ri 2 $2 

(26) 

* 2 (4 + #W + ay) + (r + a){9 2 ~ 

tt 2 (^ + tf|)(r + a) + 12 (rf + arj) 

Here and hereafter, we let the dot above a symbol denote the derivative with respect 
to r and omit the argument r of the functions rj, rj, r]', g 2 , $ for brevity. Moreover, we 
can replace g 2 with its & equivalent (see appendix). We can calculate this solution in 
the variables (x,y) if we supplement relations (13), (14), (16) with additional relations to 
obtain the complete basis of the transition between the 'r-' and '^-representations'. We 
can do this using the expressions for the quantities $2 and m 

^= 2 K', ti 2 = -$ 2 K'E', r] = K'E' — x' 2 . 

IT IT 3 

Simplifying and replacing a. i-)- ia, we obtain the compact solution 

x(aK' + K) 
V ~ aE'-E + K' 
We can analogously simplify the two other limit cases: 

2(aK' + K)(E' - K') - tt _ _ tt d 2(r + ia)^ 3 + ^3 

V ~ 2{aK' + K){E' - xK') - tt X P[At) ~ 2idr rf¥ z ' 
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_ 2(aK' + K)(E' -xK')-tt _ ir d 2(r + ia)i? 4 + ^4 

V ~ 2(aK> + K)(E> - K>) - ir P(z|r) ~ 2idr~ ■ 

We see that these solutions are mutually related by transformations of the form y i— > |; 

they therefore cannot originate from the modular transformations r i— y r + 1 and r4--. 
But we can easily see that they reflect discrete symmetries of equation (8) of the form z i— > 
{z + 1, z + t}, and this entails the above transformation for y, which is easily recognizable 
when written in the standard p(z\u, u/)-representation: 

(e-e'Xe-e") _ x 

p{z + u)-e = ± '- ^ y=-- 

p{z) -e y 

Of course, we can apply this symmetry to any Picard-Hitchin solution. We also note 
that we need the pair of functions (E, E') not only as a formal completion of the set of 
functions (K, K'): both the Hitchin and the Chazy solutions obtained in [16] and also any 
T], ^-function solutions to equation (1) are expressed in terms of (E,E'). 

4.3. The T -functions. We also presented examples of the Painleve c^Ln-forms for the 
degenerate Hitchin solutions because the T -functions are important for the theory of 
Painleve equations, being the generators of their Hamiltonians and other objects [14, 7, 19, 
18]. But the T -functions are primarily important not for constructing Hamiltonians, which 
are known for all Painleve equations, but for representing solutions in form (6) see [18] 
and the explanations in [7, pp. 740-1] in regard to constructing T -functions corresponding 
to the positive and negative residue series for solutions of the equation Vi- 

We note that other known examples (see the citations at the end of sect. 2.1) require 
the existence of the explicit transition 2-F1 1— > t because the 2-^1-series are essentially 
nonholomorphic functions with branching. Reducing them to the 'genuine' (r)-form is a 
very complicated problem, which is equivalent to inverting the integral ratios for equations 
of the Fuchsian type. In this respect, formulae are always 'final' in the uniformizing 
r-representation (sect. 5), while the transition 2-^1 > t is nontrivial and may contain 
several hypergeometric series. For instance, in the above examples, besides the integrals 
K, K' , E, E' , i.e. functions of the type 2-^1 (±5, \] we encounter the function r/ 2 (r). 
It can be shown [4, p. 263] that the hypergeometric equation related to this function is 
just the equation for the function 2-^1 (xg, §1^)) i- e - the- equations that determine the 
classic Klein modular invariant J = (g%/ 'A)(r) [1, §10], [8, §14.6.2]: 

By this means we currently have only the Hitchin class in which all the results including 
degenerations can be written explicitly. It remains to note that the found Painleve form 
admits a natural interpretation, analogous to formulae of the straight-line sections of the 
two-dimensional Jacobians in soliton equation theory. Function (21) is the canonically 
normalized meromorphic elliptic integral I(z — B\t), which has only the first-order pole at 
the point z = B and is considered on the straight-line section z — At + B = of the space 
{C x H + }, obtained as the direct product of the elliptic curve (the parameter z) and the 
one- dimensional moduli space of elliptic curves (the parameter r). 
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The T -function constructed in [15], corresponds only to Picard poles (19) and coincides 
with the function 



T S (.,,A.B) - ] {A§^ ) +B 



KWxy\ 



up to a fixed nonmeromorphic factor. Hence, there exists a second Hamiltonian corre- 
sponding to the second T -function and poles with the opposite signs of residues. This 
function is given by the formula 

T 2 (x; A, B) = TiO; A, B) A Ln { Tl(x; A> B )^ AB ) . 

We choose such a representation because meromorphic Abelian integrals (integrals of the 
second kind) of the type of functions (21) can be represented as derivatives of logarithmic 
integrals (integrals of the third kind) with respect to the parameter determining the loca- 
tion of one of the two logarithmic poles (with respect to the constant B). This may entail 
further generalizations, but we do not discuss this here. 

5. The equation V 6 and uniformization of curves 

5.1. Algebraic solutions. Having the general integral leads to several corollaries, for 
instance, to the limit cases with respect to the parameters (A,B), for example, 

„ n m , , 

At + B = -t + - (27) 

We thus obtain an infinite series of algebraic solutions of the form P(x, y) = 0, where 
P is a polynomial in its arguments. Hitchin also presented particular solutions of the 
algebraic type, and the origin of the whole algebraic Picard series was revealed in [16] 
where asymptotic regimes and monodromies related to this and other solutions and the 
character of their branching were also considered. 

It is less known that some of these results and also the very scheme for constructing 
such solutions were already contained in Fuchs' paper [24] with the same title as his 
well-known paper 2 of 1907. The procedure for constructing solutions was missed in [10] 
but the method for obtaining these solutions is obviously the same: the theorems of 
multiplication and addition for elliptic functions can be used. Moreover, Fuchs pointed 
out [ ] how to use the known Kiepert determinant multiplication formula for p{Nz) as 
a rational function of p(z) [ , p. 19], [ , p. 332]; he also wrote the Puiseux series for 
these solutions [10, §IV] and considered the monodromy of the Fuchs equations in great 
detail. Such determinant schemes can be interesting in themselves because the list of 
algebraic solutions of the Painleve equations steadily increases. We restrict ourself to the 
statement (without proof) that we can obtain effective recursion schemes for constructing 
such solutions. We do not present examples of these recursions because another immediate 
corollary from what was said above seems more interesting to us. 

From the explanations in sect. 1.1 and from explicit form (18) of the z(r)-dependence 
for the solution to equation (2), we obtain infinite series of algebraic curves P(u,v) = 



2 In a one-line footnote in [22, p. 300] Picard also claimed that the infinite series of algebraic solutions 
appears in the case where the numbers (A, B) are commensurable. 
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of nontrivial genera g > 1 with the explicit parameterizations by the single-valued 
analytic functions u(r), v(r). Because the Kiepert formula above contains a redundant 
set of functions (p, p', p", . . .), it is more convenient to pass to the language of theta 
functions, which also provides a more compact answers. In this case, it is desirable to have 
multiplication theorems directly for the functions 6^, 6[. To the best of our knowledge, 
such formulae are absent in the literature, and we therefore present them in the appendix 
(Lemma 2). For the Hitchin solutions, we can use transformation (20), but we lack the 
efficiency when acting with them on algebraic functions whose complexity rapidly increases 
as N increases. It is therefore more convenient to use the abovementioned theorems for 
the 6, ^'-functions directly in formulae (12), (18) (also see Proposition 2 below). It is 
obvious how we can apply these theorems to the solutions under investigation, and we 
therefore pass to demonstrating specific examples and their corollaries. 

5.2. Examples. Setting At + B = g, we reconstruct the Dubrovin solution [16, (4.5)] 

y 4 -(6y-4)yx + (4y-3)x 2 = 0. (28) 

Formula (3) then results in the curve of genus g = 3 of the form 

4 ( . 2 + = -3, = - = ffpy} ( 29 ) 

if we use the substitution (x = u 4 , y = —u 2 v). This curve does not belong to the series 
of classic modular equations [8, §14.6.4] nor to their closest variations [ , §7]. In turn, it 
is remarkable that because of the birational transformation 

( x 4 -y-l ( x = - ^ 



u 



Ax 
x 4 - y + 1 



2(v? + v) 



_ (v 2 + 3) 2 f 4u 2 v 

V - W 2(„2_ 1 )\ v ^2^ + 



2x 2 

equation (29) is isomorphic to the classic hyperelliptic Schwarz curve 

y 2 = « 8 + 14a; 4 + l. (30) 

In Schwarz's Gesammelte Abhandlungen [ :] this curve appears in many contexts in both 
volumes 3 , but no variants of its parametrization were known. We present the formula for 
x(t), because it admits the compact simplification 

For this curve and other curves and functions introduced by this method, we can write 
the remaining uniformization attributes (which is a purely technical problem) : differential 
equations for function y(r), the 2-^i-solutions of the corresponding Fuchsian equations, 
the representations of the monodromy group of these equations, the Poincare polygons 



3 For instance, it appears in the 1867 paper [24, I: p. 13], where the famous Schwarz derivative tf'^u) 
was introduced. 
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8], the Abelian integrals as functions of r, their reductions, if any, to elliptic integrals, 
etc. 

If we consider analogous solution (9), with the same parameters At + B = |, then we 
obtain the curve that is the Hitchin solution 

y 4 - 2 + 1) (y 2 + x)y + 6xy 2 + x 3 -x 2 + x = 0. (32) 

But this curve with the substitution x = u 4 must become isomorphic to the previously 
described curves because the initial Picard (28) and Hitchin (32) solutions have the genus 
5 = 0. In turn, formulae (3) and (9) provide an equivalent but not obvious parameteriza- 
tion of (32) 

which follows from the nontrivial $, 6, ^'-function relations' 1 . The general statement about 
the algebraic Hitchin solutions follows from formulae (3) and (10). 

Proposition 2. The algebraic Hitchin solutions admit the parameterization 

x = ^ y = ^e!\W^i^- d2 \' (33) 

where the 9 , & -functions are taken with the arguments in (27). Functions (33) satisfy 
autonomous third order differential equations of the forms 



x 3x 2 1 x 2 - x + 1 V 3 y 2 



Q(x,y) 



x 3 2i 4 2x 2 (x-l) 2 ' y 3 2 y 4 

with the explicitly calculable rational function Q(x,y) for all (n,m,N). 

Because of this proposition and the simplicity of formulae (33) , there is no need to use 
the equations P(x, y) = because the uniformization allows investigating solutions with 
arbitrarily large (n,m,N). Here, we include differentiation, the Laurent-Puiseux series, 
calculation of monodromies, the branching, the genus, plotting of graphs, etc. 

The zero genus is not a general property of the solutions under study. For example, 
under At + B = g the Picard solution has the form of a cumbersome but elliptic curve 

P(x,y)=0: 

y 12 - 50xj/ 10 + 140(x 2 + x)y 9 - 5(32x 3 + 89x 2 + 32x)y 8 

+ 16(4x 4 + 35x 3 + 35x 2 + 4x)y 7 - 60(4x 4 + 13x 3 + 4x 2 )y 6 + 360(x 4 + x 3 )y 5 (34) 
- 105xV - 80(x 5 + x 4 )y 3 + 2(8x 6 + 47x 5 + 8x A )y 2 - 20(x 6 + x 5 )y + 5x 6 = 

with the Klein modular invariant J = |||. The Weierstrass p-parameterization for (34) is 
very complicated and less adjustable for this curve, while for any r £ H + we have 

_ Kir) g|(r)g(ibj 



4 In particular, 8'i(pT + q\r) is algebraically related to the functions 8(pT + q\r), for rational p, q while 
this relation can be only differential in the generic case (see Lemma 1 in the appendix). 
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The curve genus increases rapidly. For N = 4,5, 6, 7, 8, 9, . . . the genera of the respective 
Picard-Hitchin solutions are g = 0,1, 1, 4, 5, 7, . . .. We mention here that the genus de- 
pends on the choice of the representing equation with the Painleve property. For instance, 
the same solution (34), but in Picard representation (23-24) has a nonhyperelliptic genus 
g = 3, and a simple solution of equation (24), which is birationally isomorphic to the 
lemniscate w 2 = z A — 1 (the calculation is nontrivial), already appears at At -+- B = \: 

x A u 16 - 20x 3 u 12 + 32x 2 (x + l)u 10 - 2x(8x 2 + 29x + 8)u s 

(35) 

+ 32x(x + l)u 6 - 20xu* + 1 = 0. 

Remark 2. In [ ], we advanced an assumption about a possible relation between the 
theta constants Q{u(t)\t) of general form and the uniformizing Fuchs equations. A par- 
ticular case of such a constant is just the Dedekind function 

2Ei T 

f](j) = — ie 3 6\{t\3t), 

which appears both in our formulae and in some known but isolated parameterizations of 
modular equations (mostly of the genus zero; see the bibliography in [4]). It follows from 
the above results that the Painleve equation admits an infinite class of generalizations 
in this direction using the fundamental object 0[^t + 77 |t), together with the explicit 
formulae for curves that are not necessarily solutions of the equation Ve, for functions 
on these curves, for the differential equations on these functions corresponding to the 
subgroups of r(2), and so on. 

We note the following computational fact. The Painleve curves are specific and not 
always easily accessible by various algorithms. Curve (35) is a good example. But in 
the framework of what was said above, convenient equations appear, which in particular 
include an important class of hyperelliptic curves. Even less obvious is that function (31) 
in fact corresponds to a larger curve 

z 2 = (x 8 + Ux 4 + l)(x 5 - x) . (36) 

Namely, function (31) also parameterizes this equation in the sense that z(t), determined 
from Eqs. (31) and (36) is an analytic function that is uniquely defined in its whole 
domain, i.e., in H + . This example provides one more realization of a rather nontrivial 
tower of embedded hyperelliptic curves with the universal uniformizing function X(r) (see 
[4] about this function). 

We omit both the derivation of formula (36) and other theta-functional consequences of 
what was said above. Because of their large number, they deserve a separate investigation 
that is beyond the scope of this paper. This material is presented in work [6]. 

6. Concluding remarks 

Using the analytic results in both x-representation (15), (17), and r-representation (12), 
(18) we can easily obtain several useful formulae, for example, simple expressions for the 
conformal factor F of the Tod-Hitchin metric [3, 13, 25]. This factor was expressed in 
terms of the derivatives of the theta functions with respect to the integration constants in 
[3, theorem 5.1], and the solution y(x) was expressed in terms of the double derivatives 
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of the theta functions with respect to both its arguments. By Lemma 1 the final answers 
are actually expressed in the basis i9, 0[, 0k(Ar + B\r). Another example is third-order 
differential equations for the T -functions and for the entire holomorphic functions whose 
ratio gives the solution y(x). On numerous occasions, Painleve himself indicated the 
existence and importance of these equations (see, e.g., [20, p. 1114]). They follow from 
the third-order equations for a(ar + b\r), which follow from Lemma 3 in the appendix. 

Recalling the remark in sect 1.1, we note that the known solutions and local branch- 
ings [11, 16, 13, 12] are compatible with the uniformization by the Painleve substitution, 
i.e., they generate the r-representations of the solutions in terms of the single- valued 
functions y(r). The necessary data for establishing this property for all the solutions is 
presumably contained in Theorem 2 in [18] and in results of work by Guzzetti [12]. But 
in the uniformization context, it is desirable to have the direct statement about the global 
single- valuedness of the function p[z(t)\t) for all values of (a, (3, 7, 5). This would provide 
a new important interpretation of the equation V§. Equation (1) would differ from the 
lower Painleve equations in that their solutions are globally meromorphic on the plane C; 
for the equation Vq, they are globally meromorphic on the upper half-plane IH + with the 
constant negative Lobachevskii curvature, which is the universal covering of C\{0, l,oo}. 
A broad class of algebraic solutions is then provided by automorphic functions on the Rie- 
mann surfaces with three punctures (orbifolds), while the transcendental solutions of the 
Picard-Hitchin type themselves generate another, physically meaningful class of single- 
valued functions. In this context, solutions to equation (2) do not follow automatically 
from those of (1), and integrating (2) is therefore an independent problem. 



7. Appendix: rules for differentiating the theta-function 
7.1. Jacobi's ^-functions. These functions are used in the following definitions [1]: 

e l{z \ T) = _ ie 4™ ^ ( _ 1) fc e (fc»+fc)«iT e ( 2 *+i) W i* j e ^ T) 

k=— 00 

1 • °° 

6 2 (z\t) = e 3 ™ Y, e^ +k ^ iT e^ k+1 ^ iz , 9 4 (z\r) 

k=— 00 

We also introduce the fifth independent object 

1 ■ °° 

6'Mt) =vre 3mr £ (-l) fc (2k + 1)e (* 2 +^T e (2k+iy lz 

fc= — 00 

and the i?-constants are defined as "d^ := , &k{ T ) = 0fc(O| r ) an d i?i = 0. 

Lemma 1. The five J acobi functions 01,02,03, 0a 0\ '■= -q~ satisfy the closed ordi- 
nary differential equations with respect to the variables (z, r) over the field of coefficients 



J- e fc2 --e 

k= — 00 
00 

E(-» 



2kniz 



k gfc 2 7riT ^2kmz 
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7/(t) and $ z 



cm 

dz 

mi 

dz 
dr 



d6[ 









01 




el 




0x 



7T 



(37) 



12 



4tt 9\ 



^vi + -(4 + n)\-0i 

{ 



nt 0V0, 



+ 



7T1 

T 



4#l 



+ 



7T 



7T 



■0k 



-i ^3 

4^~eJ + 7 



7T 



2+^ + ^1 + ^ 



7T 



02030 A 

~0T 



where 



8k -28 



10k - 28 



ft 



k = 1,2,3,4. 



3k- 10' ^ 3&-8 ' 
These differentiations are not completely symmetric, and the equations themselves mean 
that not only elliptic functions (the ratios 0j/0k) have the closed differential calculus 
with respect to both the variables (z, r) independently but also the same holds for their 
'constituents', the theta functions including the function e[. The details and, not less 
important, the analysis of the integrability condition for the above equations have been 
detailed in work [ ]. The basis #1,2,3,4 can be obviously closed by adding any of the 
functions 9' k , not necessarily e[. It is essential that the known equation Atti9 t = ZZ must 
be treated as a corollary of the above equations and not vice versa because the Lemma 1 
defines ordinary differential equations, while this equation is an equation of the heat kernel 
type in partial derivatives; the theta functions are just particular solutions of it. 

Lemma 2. Let n be an integer and n 1 := n — 1. Then the Jacobi functions satisfy the 
recursive multiplication formulae 



9i{2z 

9\(nz 

9 2 (nz 
9 3 (nz 
9 4 (nz 



20i (a: 



9 2 (z)9 3 (z)9 4 (z) 



i(n 1 z)9l(z)-9i(n 1 z)9l(z) 



(38) 



0| 


■9i{{n 




2)z) 






9 2 3 (n lZ )6 


!(*) - 




(n x z) 


0\ 




n 


■9 2 {{n 




2)z) 






9\{n x z)b 


!(*) + 


n 


( n i z ) 


01 






■9 3 {(n 




2)z) 






0i(n lZ )6 


S(*) ~ 




(n lZ ) 


01 





$l-9 4 ((n-2)z) 
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We obtain the multiplication formula for the function 9[{nz) by taking the derivative in 
(38), and subsequently using formulae (37). The lemma statement also holds for an arbi- 
trary complex n. 

It is natural here to assume that there are general nonrecursive expressions for these 
formulae in form of certain determinants. 

7.2. The Weierstrass functions. The modular Weierstrass functions g 2 (T), 9s( T ) an< ^ 
rj(r) are defined by the series 



T](T 



2ir z 



- oo 

24 ~ 2 



,2Ai7riT 



k=l V 1 e ) 
I k=l 



2 ( ' 



^2kir'iT 



2knir 



7T 

24 



(39) 



93 i T ) 




2knir 



7T 

432 



Lemma 3. TTie Weierstrass functions o~,£, p, p'(z\r) satisfy the nonautonomous dy- 
namical system with the coefficients r](r), g 2 (t) and wrai/i i/ie parameter z: 

da 

- _<!,,-(- 4- •>„(-•• - n - 

12- 

1 

6" 

2 



dp 
9r 



l -{p-Q 2 + 2n(zQ-l)--g 2 z 2 }<j 
l -[p' + 2{Q- zn)p + 2nQ- p 2 zj 
1 -{2{C - zr,)p' + ^{p-r,)p- 2 -g 2 } 



(40) 



--{6( p -n)p> + (C-z V )(l2p 2 -g 2 )} 

7T J 



where we use the right brace additionally to denote the differential closedness of the func- 
tions (,p,p'. Using the replacements £ h- )■ £ — 277 and <r 1— >• a exp{ — \r\z 2 } we eliminate 
the parameter z from Eqs. (40), which is equivalent to passing to the functions 9' and 9. 

Because p 1 can be expressed algebraically in terms of p, the 'essential' part of system 
(40) is its second and third equations. In turn, they are equivalent to the equations for the 
functions (,(At + B\t) and p{Ar + B\r), if we take the known relation d z C,{z\r) = —p(z\r) 
into account. Therefore, if we regard the parameters (a, 0, 7, 8) as the moduli of the 
equation Ve, then the Picard-Hitchin-Okamoto class (see Remark 1) has common moduli, 
and all the equations V& in this class are then equivalent to a single representative, which 
is system (40). 

Corollary 3. The canonical representative of the equations in the class of the 
Picard-Hitchin-Okamoto parameters is the system of two equations for the functions C,, p: 

g_i { „ +2(p+1 , )c} , £ 



~{2Cp' + 4(p-f7)p-!<7 2 } (41) 
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and its common integral 

C = ((At + B\t) - A V '(t) - Bt){t) , p = p(Ar + B\t) . 
Here p' = y / 4f? 3 — g 2 p — g 3 , and Eqs. (41) are supplemented by their corollary: 

^ = -i{6(p- 7? )p' + (12p 2 - 52 )C}. 

It is clear that we can keep either of the quantities p, p', and solutions and their 
derivatives are then rational functions of (£, p, p'). Such derivatives can then be treated 
as the p-form of Okamoto transformations (20). For instance, for the Hitchin solution 
H = p(z\t), we have 

1 P' ^ m 4H 3 -g 2 H-g 3 o tt2 a tt 1 

H = p+-*-, ^— = / 2 y3 +2H 2 + 4 V H+-g 2 . 

2 Q i dr p — H 6 

Choosing a convenient basis composed of (£, p, p') based on governing equations (41), 
we can go further and develop the inverse, i.e., the integral Painleve calculus in both 
the uniformizing (r)-half-plane and the (x)-plane. There are numerous examples, and in 
addition to the first obvious example following from the T -form of (15) 

' ydX -2Ln 





£) + 2TrA-e 1 (A£ T + B\ 






x K' ■ 9t(A§ + B\ 





x(x — 1) 

we present another, more remarkable example: 

T 

^ j{p ~ C)dr = hn9i(\AT + \B\t) - Ln7y(r) + jA 2 r. 

This example, which generates all other integrals, can be logically considered the Painleve 
(r)-analogue of the well-known integral Weierstrass relation between the meromorphic 
objects C) P an d the entire function a. Again imposing condition (27), we obtain various 
theta constants #(u(t)|t) , 9' (u(t)\t^ of the general form in the left-hand side and the value 
of the integral calculated on these constants, itself a theta constant, in the right-hand side. 
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